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' Abstract 

The V-mode polarization of the Cosmic Microwave Background is discussed in a weakly 
magnetized plasma. The VV and VT angular power spectra are computed for adiabatic ini- 
tial conditions of the Einstein-Boltzmann hierarchy. Depending upon the frequency channel 
and upon the magnetic field intensity, the VT power spectra of the circular polarization can 
^ ■ even be seven orders of magnitude larger than a putative B-mode polarization stemming 
from the lensing of the primary anisotropics. Specific programs aimed at the direct detec- 
tion of the V-mode polarization of the Cosmic Microwave Background could provide a new 
observational tool for the scrutiny of pre-decoupling physics. 
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The spectral energy density (per logarithmic interval of frequency) of the Cosmic Mi- 
crowave Background (CMB in what follows) is maximal, today, for photon energies E™ ~ 
9.2 x 10~ 4 eV whose associated wavelength A 7 ~ is of the order of 0.13 mm. Ac- 

cording to the WMAP 5-yr data the redshift of hydrogen recombination can be estimated 
approximately as z ICC = 1090.51 ± 0.95 [TJ corresponding to the conformal time§ r rec ; at 
r rcc the maximum of the CMB is z rec E^ ~ eV corresponding to a physical wavelength 
A 7 (T rcc ) = A 7 /z rcc — 0.12 fim. Prior to recombination the electrons and the ions have kinetic 
temperatures which are comparable with the temperature of the photons, i.e. (1 + z rec )T j0 
(where T 70 = 2.725 K). The small difference between electron and proton temperatures is 
controlled by the ratio between the Hubble rate H and the Coulomb rate Tcoui which is 
(9(10~ n ) at r rec . The global charge neutrality of the plasma combined with the baryon 
asymmetry implies that the electron and proton concentrations are equal and both of the 
order of 10~ 10 n 7 where n 7 is the comoving photon concentration. Prior to recombination 
the plasma is cold: the electron and proton masses are both much larger than the kinetic 
temperature of the corresponding species. Consider the physical situation when, prior to 
recombination, the plasma is supplemented by a magnetic field whose typical inhomogeneity 
scale is at least comparable (and possibly even larger) than the Hubble radius H^ 1 at the 
corresponding epoch. Since the wavelengths of the scattered photons are minute in com- 
parison with the Hubble radius (i.e. A 7 (r rcc ) <C H^ 1 ) the magnetic field gradients can be 
ignored, in the first approximation, when computing the photon-electron (and photon-ion) 
scattering. The gradient expansion on the magnetic field strength was termed long ago by 
Alfven guiding centre approximation [2]. 

Having introduced t\ and e 2 as two mutually orthogonal directions (both perpendicular 
to the direction of propagation of the radiation), and recalling the standard definitions of 
the Stokes parameters [3] it can be easily shown that 

I = \E-e 1 \ 2 + \E-e 2 \ 2 , V = 2 Im[{E ■ e x )*{E ■ e 2 )], (1) 

are both invariant for a rotation of i\ and e 2 on the plane orthogonal to the direction of 
propagation of the radiation. For the same two-dimensional rotation, (Q ± ill) transform as 
a function of spin weight =p2 on the two-sphere [3]; this observation leads, after some algebra, 
to the known form of the E-mode and B-mode polarization jSJ. If a large-scale magnetic field 
is present and if, concurrently, the spatial curvature does fluctuate over large scales, then 
the power spectra associated with the brightness perturbations of V will not be vanishing 
and shall be defined, in what follows, V-mode power spectra in analogy with the B-mode 
and E-mode power spectra characterizing the linear polarizations. 

2 The conformal time coordinate r will be used throughout; in terms of r the background metric g^ u will 
be chosen as conformally flat i.e. g — a 2 ^)^^ where rj^ — diag(l, —1, —1, —1) is the Minkowski metric. 
The current observational evidence pQ suggests, indeed, that the spatial curvature can be neglected at the 
recombination epoch, at least in the framework of the concordance model, i.e. the ACDM paradigm where 
A stands for the dark energy component and CDM stands for the cold dark matter contribution. 
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To compute the induced V-mode polarization the evolution of the brightness pertur- 
bations must be written in the case when the photons scatter electrons in a magnetized 
background. In the elastic e-7 scattering occurring in a cold plasma the recoil energy of 
the electron can be neglected [B]; photons impinging on electrons and ions in a weakly mag- 
netized medium can be described, as usual, in terms of a scattering matrix connecting the 
outgoing to the ingoing Stokes parameters (see, e.g. [3j). The latter step will lead, after 
angular integration, to the evolution of the various brightness perturbations^ 



A; + ifc//(Ai + 4>) + e'Ai = xp' + e'Ci(u, k, //), 
Ap + ik/jAp + e'Ap = e'C P (uj, k, 
Ay + ikfxAy + e'A v = e'Cy(u, k, //), 



(2) 
(3) 
(4) 



where the prime denotes a derivation with respect to the conformal time coordinate r while 
4> = 8^ goo/ '(2a 2 ) and ipSij = 5^gij/(2a 2 ) are the scalar fluctuations of the metric whose 
relation to the curvature fluctuations can be expressed, in the longitudinal gauge, as 



n = -tp- 
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note that the relation of % to the Hubble parameter is simply given by aH = TL; the source 
functions appearing in Eqs. 02]), flH]) and (j3J) are given by 
1 
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where the dependence upon the frequency of the observational channel arises through the 
functions 
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3 In what follows e' — ax c nga 7C is the differential optical depth, no = no/a 3 is the physical concentration 
and <7 7C = (8/3)7r(e 2 /m c ) 2 . The quantity u is simply the projection of the Fourier wavevector on the direction 
of the photon momentum. 
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The plasma and Larmor frequencies for electrons and ions are denoted, respectively, by 
ubc) and (co v u w b0- CMB experiments operate for angular frequencies which are 



to 



typically larger than the Larmoijj and plasma frequencies of the electrons at recombination; 
it is therefore legitimate to expand the source functions of Eqs. (j2]), ((3]) and (00) in powers 
of f e (u)) = (wbc/w) as well as in powers of (m e /m p ); the result of this double expansion can 
be written as 
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where Pe(p) is the Legendre polynomial of £-th order and where Sp(k, r) = (A P0 + A P2 + A I2 ) 
is the standard source term for the E-mode polarization when f e (u>) = 0. Indeed, in the limit 
f e (u>) — > 0, Eqs. (EJ), (EJ) and (jlj) reproduce the standard results for the evolution equations 
of the scalar brightness perturbations. The source functions obtained in Eqs. (|T0|) - (|T2|) 
are derived by integrating the angular dependence of the ingoing Stokes parameters in full 
analogy with what happens in the case when the magnetic field is absent [3] (see also [7] for 
further details as well as [8] for slightly different perspectives on magnetized photon-electron 
scattering). Note that 
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where UBe is the Larmor frequency and B n denotes the uniform component of the comoving 
magnetic field intensity which is treated within the guiding centre approximation (see [9] for 
the description of magnetized plasma prior to recombination). The numerical solution of the 
system is greatly helped by exploiting systematically the integration along the line of sight 
[5] for all the brightness perturbations. From Eq. (jlj) 
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4 Note that u; pc /w p i = ^m p /m e where m p is the proton mass; ui = I'kv denotes throughout the angular 
frequency. 



where /C(r) is the visibility function and were (tq — r) is effectively the (comoving) angular 
diameter distance in a spatially flat geometry. To zeroth order in the tight-coupling expansion 
Eq. (|T4|) allows to evaluate the V-mode polarization, i.e. 

A v (A;,/i,7b) = I r K{T)f c {ujy- lllk{T0 - T) PM Ai (fc,r)dr. (15) 

O JO 

where Ato(A;,t) is the monopole of the intensity computed to lowest order in the tight 
coupling expansion, i.e. when the baryon velocity coincides with the dipole of the intensity 
of the radiation field. To lowest order in the tight coupling approximation the CMB is 
circularly polarized provided a large-scale magnetic field is present. The linear polarization 
is generated to first-order in the tight-coupling expansion but is larger than the V-mode 
polarization because of the smallness of / e (u;). More details on this semi-analytic discussion 
can be found in [7] . The evolution of the monopole of the intensity Ai can either be studied 
in the tight-coupling limit or it can be solved numerically. In both cases the initial conditions 
will be chosen to be adiabatic % 

Aio(^,r rcc ) = ^±^ TZ,(k), Mk) = (l + \Rv)Mk), (16) 

where R v is the fractional contribution of the massless neutrinos to the radiation background 
and H*(k) denotes the curvature perturbations prior to equality and for typical scales larger 
than the Hubble radius at the corresponding epoch. Since prior to equality H = 1/t we 
shall also have, from Eq. that TZ* (k) = —ip*(k) — <f)*(k)/2. For large angular scales 
(i.e. i < 50) the visibility function can be considered to be sharply peaked at recombination 
and it is in practice a Dirac delta function. The V-mode autocorrelation (i.e., for short, VV 
power spectrum) and the cross-corrrelation between polarization and temperature (i.e., for 
short, VT power spectrum) can then be computed analytically in this regime and the result 
|f| can be written as: 

(17) 



75 (4^ + 15) 2 JeK ' "-\k 



5 See, e.g., [§]. The present considerations can be straightforwardly generalized to the case of non-adiabatic 
initial conditions |10j . For a discussion on the peculiar features of the adiabatic and non-adiabatic initial 
conditions of the Einstein-Boltzmann hiererchy see, for instance, [IT] . 

6 For simplicity the angular power spectrum shall be denoted as g1 X = £(£ + l)Cl XY ^/(27r). The VV 
and VT power spectra are the analog of the EE and TE power spectra arising in the case of the linear 
polarization. The temperature (related to the I Stokes parameter) and the circular polarization (related to 
the V- Stokes parameter) are both invariant under a rotation orthogonal to the direction of propagation of 
the radiation, as stressed after Eq. (JTJ . It is therefore natural, in a first approach to the problem, to consider 
the TT, VT and the VV power spectra. Furthermore, the VT correlations are larger than the the VE and 
VB correlations. 
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where k = (r — r rcc ) 1 and k p = 0.002 Mpc 1 is the pivot scale while An is the amplitude 
of the power spectrum of curvature perturbations at k p (in the concordance paradigm and 
in the light of the WMAP data alone An = (2.41 ± 0.11) x 10~ 9 ); the functions Z\ 



<w), 



and zj VT ^(n s ) appearing in Eqs. ( TT7|) and (|TS|) are nothing but: 

+ \)\U{t + 1) - (n s - 1)K - 2)(n s - 4)]r(3 - n B )Y(t - f + f ) 



xf v) K) 
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^+l)(2-n s )r(2-f)r(^ + f -1 
A \n s ) = 7 \ 7 7 V , ^, (19) 



and arise as analytically solvable integrals of products of spherical Bessel functions and of 
their derivatives. Before discussing the relevant numerical results over small angular scales 
it is appropriate to mention here that circular polarization is often invoked as the result 
of the Faraday conversion^ of linearly polarized radiation [12]. For the latter mechanism 
to operate, relativistic electrons must be present in the system and this can happen only 
as a secondary effect when CMB photons pass through magnetized clusters (see [12], last 
reference); this is however not the idea pursued here since the pre-decoupling plasma is cold 
and electrons are deeply non-relativistic. The V-mode polarization, as we showed, is induced 
by the magnetized plasma itself thanks to the presence of (adiabatic) curvature perturbations 
in the system. Absent one of these two components the VT and VV power spectra would 
vanish. If the initial conditions would not be adiabatic the V-mode polarization would still 
be present but with different physical features which will depend upon the specific non- 
adiabatic solution [TTJJ E]. For smaller angular scales (i.e. £ > 100) it is mandatory to 
integrate numerically the system across decoupling. Some of the results are summarized in 
Fig. [TJ The thin lines in both plots denote the V-mode autocorrelations while the thick lines 
denote the cross-correlation of the circular polarization anisotropies with the temperature 
inhomogeneities. The signal is larger for low multipoles and its shape reminds a bit of the 
temperature autocorrelations induced by the tensor modes of the geometry which reach their 
largest value for small i and decay exponentially for i > 90. Defining as rx the tensor to 
scalar ratio at the pivot scale k p (TJ |TT] , for rx = 1 the TT correlations induced by the tensor 
modes would be O(10 3 ) (/iK) 2 while the VT correlations are O(10~ 5 ) (/uK) 2 for the choice of 
parameters of Fig. [T] (see, for instance, plot at the right, dot-dashed curve). For the same 
choice of parameters the VT power spectra are of the order of the B-mode autocorrelation 



7 Faraday conversion (typical of relativistic jets) should not be confused with Faraday rotation. In the 
presence of relativistic electrons linearly polarized radiation can be Faraday converted into circularly polar- 
ized radiation [T2] . Faraday rotation is instead a rotation of the polarization plane of the (linearly polarized) 
radiation: in practice it can convert E-modes into B-modes but it does not lead to circularly polarized pho- 
tons (see, e.g. third and last references of [S])- Faraday conversion and Faraday rotation have also different 
dependences upon the magnetic field intensity and upon the frequency. 
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Figure V. In the plot at the left the VT and the VV angular power spectra are reported for 
a fixed value of B u (i.e. 1 nG) but for different values of the comoving frequency. In the 
plot at the right the comoving frequency is fixed to 10 GHz but the magnetic field strength 
increases. The thin lines denote the VV correlations while the thick lines denote the VT 
correlations. In the plots, on both axes, the common logarithm of the corresponding quantity 
is reported. 

induced by the weak lensing of the primary anisotropies (i.e. (9(10~ 5 )(/iK) 2 , see also left 
plot of Fig. [2] and the discussion below). By shifting the observational frequency the VT 
correlation can be even larger [T3t [TJ]. The B-mode autocorrelation induced by the tensor 
modes of the geometry is typically larger, both than the V-mode polarization and than the 
BB spectra from lensing. For rx ~ 0.1, the BB angular power spectrum of the tensor modes 
of the geometry is (9(10~ 2 )(yuK) 2 for I ~ 90 corresponding to angular separations of roughly 
2 deg. In Fig. [2] (plot at the left) the EE power spectrum stemming from the best fit to the 
WMAP 5-yr data alone is illustrated with a dashed line and compared, in the same plot, with 
the BB angular power spectrum arising from the lensing of the primary anisotropies (thin dot- 
dashed line) as well as with the V-mode autocorrelation (full thin curve at the bottom). The 
B-mode autocorrelation stemming from the tensor modes in the case rx = 0.1 is illustrated 
with the thick line. Always in Fig. |2] (plot at the right) the TE and the VT correlations are 
compared. Both in Figs. [T] and 12] the cosmological parameters are fixed (as indicated in the 
title of each figure) to the values of the best fit stemming from the WMAP 5-yr data alone in 
the light of the concordance model. In Fig. [2] the frequency of the channel has been taken of 
the order of 30 GHz. Even if the latter frequency is already rather low, it would be desirable 
to reduce it even more and to conceive spectropolarimetric measurements in the range of the 
GHz. The challenge of detecting the CMB radiation at low frequencies is neatly described 
in Ref. [13] where a set of absolute radiometers is employed in different channels at 0.6, 0.82 
and 2.5 GHz (see also [H] for earlier results along the same theme). As specifically discussed 
also in analytic terms (see Eqs. (fTT|) - (TT8"]) ) the VV and VT power spectra are sensitive to 
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Figure 2: The V-mode power spectra are compared with the linear polarizations. As in Fig. 
[Don both axes the common logarithm of the indicated quantity is reported. 

the underlying cosmological parameters, to the initial conditions of the Einstein-Boltzmann 
hierarchy as well as to the magnetic field parameters. For illustration the concordance model 
supplemented by adiabatic initial conditions has been considered. The maximal intensity of 
the comoving magnetic field has been taken to be of the order of the nG. This is the range 
of current bounds stemming from the simultaneous analysis of the measured TE and TT 
power spectra (see [9], first and second reference). Larger magnetic fields would distort the 
acoustic oscillations in the TT power spectra. Low frequency instruments could make the 
difference for scrutinizing a potential V-mode polarization. In this respect the results and 
the techniques of p2] (as well as the earlier results of [33]) could be probably revisited in 
the light of the considerations developed here. It has been shown that the VT correlation 
for a comoving magnetic field from 5 to 10 nG can be as large as 10~ 5 (/iK) 2 at 10 GHz 
for I < 20 (i.e. large angular separations). This means that for frequencies O(MHz), the 
resulting signal could be even 6 or 7 orders of magnitude larger than a putative B-mode 
signal from gravitational lensing (see, e.g. Fig. [2j thick dashed curved in the left plot). It 
has been demonstrated that the study of circular dichroism is not more forlorn than other 
signals which are often invoked as conceptually important to consider but observationally 
difficult to assess. The systematic effects plaguing the measurements of the V-mode power 
spectra differ from the case of linear polarizations. Wether or not they are less severe depends 
also upon the features of the instrument and on the specific frequency band. The author is 
grateful to G. Sironi, M. Gervasi and A. Tartari for stimulating discussions. 
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